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ABSTRACT 

The reionization epoch of singly ionized helium (He II) is believed to start at redshifts z ^ 3.5- 
4 and be nearly complete by z ~ 2.7. We explore the post-reionization epoch with far-ultraviolet 
spectra of the bright, high-redshift quasar HS1700-I-6416 taken by the Cosmic Origins Spectrograph 
(COS) on the Hubble Space Telescope, which show strong He II (A303.78) absorption shortward of 
the QSO redshift, zqso — 2.75. We discuss these data as they probe the post-reionization history 
of He II and the local ionization environment around the quasar and transverse to the line of sight. 
We compare previous spectra taken by the Far Ultraviolet Spectroscopic Explorer to the current COS 
data, which have a substantially higher signal-to-noise ratio. The Gunn-Peterson trough recovers at 
lower redshifts, with the effective optical depth falling from Toff ~ 1.8 at z ~ 2.7 to ~ 0.7 at 
z ~ 2.3. We see an interesting excess of flux near the He II Lya break, which could be quasar line 
emission, although likely not He II Lya. We present spectra of four possible transverse-proximity 
quasars, although the UV hardness data are not of sufficient quality to say if their effects are seen 
along the HS1700 sightUne. 

Subject headings: galaxies: active — intergalactic medium — quasars: absorption lines — quasars: 
individual (HS1700-I-6416) — ultraviolet: galaxies 



1. INTRODUCTION 

The full reionization of intergalactic helium, from He II 
to He III, was likely delayed compared to that of hydro- 
gen because of the need for higher-energy photons than 
stars provi de. While hydrogen was likely rei onized at 
z > 6 (e.g.. iFan et all[2006: Shull et al.ll2012ai) . quasars 
may have created sufficien t photons with E > 4 ryd to 
reionize helium at z ~ 3-4 (jFurlanetto fc Ohll2008| ). This 
is consistent with redshift estimates from the intergalac - 
tic medium (IGM) temperature (e.g.. lBecker et al"1l2011l ). 
which increases noticeably during helium reionization, as 
well as estimates from the redshift evolution of the He II 
Gunn -Peterson optical dep th fe.g.. lSvphers et aT]|2011al 
20ia IWorseck et all 1201 it but see IDavies fc Furlanettd 
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The Gunn-Peterson test, in which a non-trivial ion 
fracti on creates a trough by li ne absorption at every red- 
shift (|Gunn fc Petersonlll965D . is the most direct test of 
the later stages of helium reionization. Unlike for hy- 
drogen, this test is useful for helium even well before 
the completi on of reionization, up to a He II fr action 
SHcii ~ 0.1 (|McQuinnl[2009tlSvpliers et al.ll2011bf ). and 
indicates that He II reionization completes at z ~ 3 or 
just below. There is noticeable sightline variance, how- 
ever, and some sightlines have x neii ^ 0.01 down to 
redshift z = 2.7 (|Shull et al.ll20Toh . 

HS1700-F6416 (henceforth HS1700) is an extremely lu- 
minous (r = 15.9) quasar at redshift zqso = 2.75. De- 
spite having several pa rtial Lyman- limit systems (pLLS; 
IVogel fc Reimerd 119950 . it is still bright in the far UV, 
with flux down to the He II Lyg br eak at 303.78 A in 
the rest frame (jPavidsen et al.l [19960 . This makes it a 
"He II quasar," one of the few percent of z ^ 3 quasars 
that have visible flux at the He II Lya break owing to 

|david.syphers@colorado.edu, micliael.sliull@colorado.edu| 



unusually low integrated H I column densities along their 
sightlines. 

In recent years, cross-matching quas ar catalogs with 
Gala xy Evolution Explorer (GALEX; iMorrissev et al.l 
l2007f ) UV imaging has lead to a vast increa se in the 
number of known He II qu asars (jSvphers et a l. 2009a, 1^, 
120121 : IWorseck et"all 120 lit ). However, He II quasars 
that are "very bright," with far-UV flux /f^v > 
10~^^ erg s~^ cm~^ A~^, remain very rare. The only 
three known, in order of brightness, are HE2347— 4342 
{Reimers et al. 1997; Shull et al. 2010), HS1700, and 
4C57.27 (jSvphers et al.ll2O120 . HS1700 is the lowest red- 
shift of these, and with no black Gunn-Peterson trough, 
it probes the helium post-reionization epoch. 

In Section [2] we present a new Hubble Space Telescope 
(HST) far-UV spectrum of HS1700, and consider instru- 
mental effects that might impact the science. In Section[3] 
we calculate the Gunn-Peterson optical depth seen in 
the HS1700 sightline, the subject of many previous ob- 
servations, and we discuss the most relevant of these in 
Section |4l Although the new observations presented are 
superior in some respects, there are resolution and wave- 
length coverage advantages in some of the older data sets, 
which allow us to better interpret the current data. We 
discuss possible quasar line emission near the He II break 
in the context of comparison to the Far Ultraviolet Spec- 
troscopic Explorer (FUSE) spectrum. He II observations, 
when combined with H I spectroscopy, uniquely allow de- 
termination of the UV background hardness throughout 
most of the IGM, and we discuss this in Section [5l In 
Section[S]we discuss proximity effects, both from HS1700 
itself and also transverse proximity effects from lower- 
redshift ionizing sources. We conclude in Section [3 and 
include a detailed discussion of our modeling of the spec- 
trum background in the Appendix. 

2. OBSERVATIONS 
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Figure 1. The COS G140L spectrum of HS1700, with a power- 
law continuum fit overplotted in blue, and continuum uncertainty 
in cyan. The COS data are binned to 0.56 A (about one resolu- 
tion element). The dashed vertical red line shows the expected 
wavelength of the He II Ly« break based on the quasar redshift, 
z = 2.748 ± 0.005. Continuum fits i ncluding partial Lyman-limit 
systems are discussed in Section 14.31 and shown in Figure |6] 

HS1700 was observed for 15,705 s on October 10, 
2011, using the Cosmic Origins Spectrograph (COS; 
iGreen et al.|[20T2t ) aboard HST. We used the G140L grat- 
ing (resolving power R ~ 2000-3000) to observe to wave- 
lengths as short as possible. The spectrum is shown in 
FigurelU including a continuum fit discussed in Section[31 
The data were taken in TIME- TAG mode, which allowed 
us to make a cut avoiding those times with evident geo- 
coronal O I A1304, resulting in 5900 s of data during or- 
bital night. In our plots and analysis, night-only data are 
used in wavelength regions contaminated by geocoronal 
emission, while all data are used elsewhere, because there 
is no evidence for a noti ceable scattered-light background 
from geocoronal lines (jSvphers et alj |201^ . The data 
are processed with CALCOS 2.18.5, using a restricted 
extraction window, pulse height amplitude cuts, and an 
improved background subtraction method. These modi- 
fications of the standard extraction method are described 
in detail in the Appendix; they noticeably affect only 
data where the background contributes significantly to 
the spectrum (A < 1100 A). 

We also require optical data to examine the H I cor- 
responding to He II absorption^ and for thi s we use the 
Keck HIRES spectrum from iSimcoe et all (|2002[ ). This 
observation gave data with a resolution of 6.6 km 
and a median S/N of 110 per resolution element. 

The redshifts quoted in the literature for HS1700 vary 
significantly. At z 3, the prominent quasar lines 
seen in optical spectra are high-ionization UV emis- 
sion (including H I Lya and C IV A1549), which can 
be shifted by many hundreds or even thousands of 
km s~^ from the system ic redshift of the quasar (e.g., 
iVanden Berk et al.ll200l . Because low-ionization metal 
lines tend to trace the systemic redshift much better, 
we measure O I A1304 and C II A1335 to find zqso = 
2 . 745ib0.004 (using the [0 II I] - corrected rest wave l engths 
of IVanden Berk et al.l [200l . iTrainor fc Steidell ([2012h 
adopted the slightly higher redshift zqso = 2.751 ±0.003 
based on Mg II A2798 in a Palomar/TripleSpec spectrum 
and H7 in a Keck/NIRSPEC spectrum. (The best lines 



for quasar redshift determination, H/3 and [O III] A5007, 
are unfortunately blocked by atmospheric absorption at 
this redshift.) The He II break is consistent with both 
of these redshifts, particularly given the wavelength cali- 
bration uncertainty we discuss in section [2Tl The differ- 
ence between the two redshifts is only 1.2a, so we adopt 
the average zqso = 2.748 ± 0.005 for this paper. The 
error is unfortunately not easily quantified, given the 
systematic uncertainties involved. Our adopted error is 
slightly larger than the statistical error in accommoda- 
tion of this. We note that small redshift variations have 
no impact on our analysis of the He II Lya forest, and 
matter only for discussion of the line-of-sight proximity 
effect fSection l6.1[) . The onset of the hydrogen Lya forest 
at z = 2.7442 is consistent with our adopted redshift, and 
clearly rules out the significantly lower zpso ~ 2.72-2.73 
used in some earlier works Ce.g.. iReimers fc Voge]|[l993l : 
iFechner et al.ll2006bl : MTei^miW . 

2.1. COS Wavelength Calibration 

The G140L/1280 central wavelength setting uses both 
segments of the COS FUV detector, with a substantial 
gap in between. Segment A covers A > 1270 A, while seg- 
ment B covers A < 1170 A. COS wavelength calibration 
is accomplished using the spectrum of an onboard PtNe 
lamp, but there are no useable lines produced on seg- 
ment B for G140L. Therefore the segment B wavelength 
solut ion is taken entirely from segment A ([Oliveira et al.l 
120101 ). We have checked the wavelength calibration of 
COS segment A using strong ISM absorption lines, which 
show offsets of ~30 km s^^, substantially smaller than a 
resolution element (the latter is ~100-150 km s~^). 

Segment B does not appear so well aligned. The only 
strong ISM lines on segment B outside the He II Lya 
forest are from Fe II, the strongest of which should be 
Fe II Al 144.94 A. This line appears to be observed at 
1142.96 ± 0.05 A; although there are other Fe II lines 
closer to this wavelength, no absorption lines are evident 
in the spectrum at longer wavelength. Another method 
of determining a wavelength offset is to cross-correlate 
the He II and H I forests, because IGM density fiuctua- 
tions should on average connect the two. A clear peak 
in the cross-correlation also suggests that the UV data 
need to be shifted to higher wavelengths, but by about 
3.4 A rather than 2.0 A. Unfortunately, the low resolu- 
tion and modest S/N of the COS data preclude a firm 
determination of the wavelength offset between optical 
and UV data. 

Cross-correlation of COS data with the independently 
calibrated FUSE spectrum also gives a strong peak at a 
2.0 A shift of the COS data, and the FUSE spectrum 
shows ISM Fc II A1143.2 A and 1144.9 A to be within 
0.1-0.2 A of their expected locations. This magnitude of 
shift also makes the He II break most consistent with the 
adopted redshift. We therefore use the COS segment B 
spectrum shifted by 2.0 A as the standard version in 
our plots and discussion, but caution that this is not an 
entirely firm determination. 

3. EVOLUTION OF HE II EFFECTIVE OPTICAL DEPTH 

In order to determine the optical depth of IGM He II, 
we need to extrapolate the quasar continuum to lower 
wavelengths. iSvphers et al- (,2012.) found that variations 
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in smoothing length and other parameters dominate the 
error in this extrapolation of COS data. We follow their 
method of varying these values and taking the resulting 
range of continua as our uncertainty. We fit a power 
law to the COS FUV data, adopting E{B -V) = 0.026 
(iSchlegel et al.lil998'l . Ry =3.1, and the extinction curve 
of iFitzpatricM (|1999f ) . The COS spectrum, with its esti- 
mated continuum, is shown in Figure[TJ Small uncertain- 
ties in the reddening are n egligible, since, for example , 
using E{B -V) = 0.022 (ISchlaflv fc FinkbeineH [2?)Tl[) 
changes our observed continuum estimation by only 2% 
at 1000 A. Lar ger uncertainties in the reddening, or 
the extinction curve are systematics we do not account 
for in our errors. 

While all we require for calculating r^ff is an extrapo- 
lation of the continuum below the He II break, HS1700 is 
one of the very few quasars for which we have continuous 
spectral coverage from the FUV to the optical. It is thus 
possible to actually measure the intrinsic power-law in- 
dex, taking into account the many pLLS. This does not 
greatly affect the extrapolated values, but is of interest 
for determining the true extreme-UV (EUV, A < 912 A) 
spectral index. The small impact on T^ff discussed in 
detail in Section 1431 

Detector background is important whenever the back- 
ground counts are significant compared to the source 
counts. This happens with faint targets, high optical 
depth regions, and — most relevant to HS1700 — at wave- 
lengths where the instrumental throughput is very low. 
We are little affected by the background uncertainty at 
A > 1100 A (He II Lya redshift z > 2.62), but at 
shorter wavelengths HST/COS has very poor through- 
put, dropping from an effective area of nearly 2000 cm^ 
at 1170 A to only 100 cm'^ at 1090 A and 10-20 cm^ at 
A < 1070 A (jMcCandliss et al.|[20lol ). The defauh CAL- 
COS background subtraction is problematic for a num- 
ber of reasons, notably the problem of "y-dip," where re- 
duced sensitivity in the primary science aperture (PSA) 
on the detector leads to a lower background in the PSA 
compared to the region where the background is mea- 
sured. For a substantially more detailed discussion of the 
COS background an d the way CALCOS handles it, see 
iSyphers et al.l (|2012t ). We have dramatically improved 
background determination for COS data, and discuss our 
method in detail in the Appendix. 

One additional concern for the calculation of t^s is 
that by necessity one usually assumes that all absorp- 
tion is due to He II. Hydrogen contamination is unlikely 
to be a major concern, as only weak higher-order lines 
of relatively rare low-redshift systems would affect this 
spectral region, but metal absorption can be noticeable. 
HS1700 is again fairly unusual among He II quasars be- 
cause it has a very well-studied high-resolution optical 
spectrum, and much effort has gone into modeling its 
metal systems. iFechner et al . (2006a) performed detailed 
modeling of the HS1700 sightlinc, predicting the metal- 
line spectrum seen in the far-UV based on observed ab- 
sorption systems and Cloudy modeling. Accounting for 
metal contamination is important when doing line fit- 
ting of the FUV spectrum (Fcchncr et al. 2006b), but 
it does not have a large effect when finding the effec- 
tive optical depth over large bins. We nonetheless do 
account for their predicted metal spectrum (Fig. 7 in 



IFechner et al.|[20"06al which also includes Galactic H2), 
as it can shift TeS by 0.1-0.15 overall, and more in some 
redshift bins. (Their model does not cover z < 2.292 or 
2.562 < z < 2.578 since it was made to cover FUSE LiF 
data only.) However, we caution that such predictions 
are not possible for many He II quasars; we therefore 
show Toff(z) in Fig. [5] with and without modeling this 
absorption. Theoretical models do not yet depend on 
such small shifts, however. 

The large-scale evolution of the He II Lya optical depth 
is shown in Figure [21 As one would expect, the aver- 
age optical depth decreases toward lower redshift. The 
COS data show a noticeably higher Tgg than was c alcu- 
lated for the FUSE data in IFechner et al.l ()2006bD : the 
COS data on HS1700 give average optical depths much 
more line with th ose seen in HE2347 at similar redshifts 
(jShuff et al.ll2010l ). Note that many previous He II stud- 
ies picked specific and irregular bins to capture what 
they felt were structures in Tes jz) (e .g., lHeap et al.ll2000l : 
iZheng et al.|[20p: iShuff et al.|[20ia ). Although theoreti- 
cal models have been coi npared to these values for want 
of any other data (e.g., iDixon fc Furlanettoll200l . we 
feel this practice is not useful in general, and bins should 
be regular and made without regard to specific sightline 
structure to enable comparison to models (which, after 
all, are trying to reproduce average behavior, not specific 
structures). However, this specifically addresses broad 
reff(z) evolution measurements. A promising method 
that should examine regions chosen a posteriori is us- 
ing "dark gaps " for helium the way it has been done for 
hydro gen (^e.g.. lPaschos fc Norman|[2005HGallerani et al.l 
I2OO6I) . 

Our uncertainties incorporate the full range of accept- 
able power-law fits to the data, but these fits do not 
include any pLLS. If we consider the best- fit contin- 
uum including pLLS, the Tcff values are slightly lower, 
as shown by the asterisks in Figure O Such a continuum 
fit would marginally strengthen the case for increasing 
Tcff at higher redshift, but in every case the points lie 
within the 68% confidence intervals for the points from 
pure power-law fits. 

In addition to broad evolution of the optical depth, 
the detailed small-scale evolution of Tcs{z) tells about 
the IGM density and local ionizing sources. We consider 
this in Section [5] 

4. PRIOR OBSERVATIONS AND POTENTIAL LINE 
EMISSION 

As one of the FUV-brightest high-redshift quasars on 
the sky, HS1700 has been observed many times before. 
It is bright enough to have been observed with the Inter- 
national Ultraviolet Explorer (lUE; Rcimers et al. 1989). 
It was observed earl y with HST using the Faint Object 
Spectrograph (FOS; iReimers fc Vogell 119931 ). which did 
not have coverage down to the He II Lya break, but 
showed that it had flux extending though the FUV. A 
Hopkins Ultraviolet Telescope (HUT) o bservation con- 
firmed that it was a He II quasar (jPavidsen et al.l 
.1996) . Subsequent observations with the Goddard High- 
Resolution Spectrograph (GHRS) and the Space Tele- 
scope Imaging Spectrograph (STIS) were not put to use 
for helium-related science, but the FUSE spectrum (cal- 
ibrated by low-resolution STIS) made it one of only two 
quasars to date with resolved He II Lya forest observa- 
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Figure 2. The broad-scale evolution of the He II Lya effective 
optical depth with redshift along the HS1700 sightline. We show 
confidence intervals on r^s including both 68% Poisson uncertainty 
on observed counts and the range of possible extrapolated continua 
shown in Fig. [T] We a c count for metal contamination using the 
model of Fe chner et al.l Il2006al) . as discussed in the text. The al- 
ternate continuum includes many pLLS, as discussed in Section [4.3l 
and shown in Figure [B] 

tions. STIS covered the FUV and NUV, thereby allowing 
a fit of the entne spectrum, and the FUSE observations 
remain the highest resolution to date, albeit at low S /N. 

To begin we make a bri ef comparison to an olde r data 
set, the FOS spectrum of IVogel fc Reimerj (| 19951 ). We 
recommend that all FUV line identifications in this paper 
be confirmed in the COS spectrum, which has twice the 
resolution and higher S/N. The FOS spectrum may be 
suspect in some areas, with systematics larger than the 
quoted uncertainties. For exampl e, the He I A584 absorp- 
tion line re ported at z = 2.1678 (|Reimers fc VogelllTOQl 
IVogel fc R eimcrs 19951) does not exist in the COS spec- 
trum, although the neighboring Galactic Al III A1854 
clearly does. Locations for claimed He I A584 at z = 
2.290 and z — 2.315 do show real lines in the COS spec- 
trum, which is important, as HS1700 is the only sightline 
with confirmed detections of this line in the IGM. 

4.1. Possibility of EUV Quasar Emission Lines 

One puzzle in EUV quasar spectra has been a lack of 
clear He II Lya emission in most quasars, despite clear 
photoionization model predictions that this line should 
be ev i dent in the broa d line region (jSvphers et al.ll2011al 
[20T1 ILawrencdlMl. The first cl aim of clearly seeing 
this line was lOavidsen" etall (IT996I ) using HUT data on 
HS1700. However, no such line was seen in FUSE data 
(Figure [3]) . It seems clear from Figures [1] and [3] that 
no reasonable continuum fit easily explains the excess of 
fiux near the He II Lya break. It is unfortunate that the 
COS spectrum has a gap just redward of this rise in fiux, 
allowing for the possibility that flux misalignment is the 
source of the problem (albeit a very unlikely one, since 
fiux matching across segments has not been problematic 
with other data). 

A fiux excess is seen above the continuum fit at 
1150 ±2 A, approximately 11 A redward of the expected 
location of He II Lya. In the rest frame, the peak occurs 
at 306.8±0.5 A (-3000 km s~i redward of He II A303.78). 
Quasar broad emission lines are often notably blueshifted 
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Figure 3. The He II Lya forest along the HS1700 sightline as 
seen in COS and FUSE. The COS data are binned to 0.32 A (just 
over half a resolution element), and the FUSE data are binned to 
0.30 A (about six resolution elements). The FUSE data are no t 
convolved with the COS line-spread function (LSF; IKrisd [20111 1: 
one can therefore see the advantage of higher resolution, in par- 
ticular being able to identify small black regions. The COS data 
show the substantial S/N improvement this detector has at longer 
wavelengths, as well as the much better characterized background. 
The dashed vertical red line shows the expected wavelength of the 
He II Lya break based on the quasar redshift. Geocoronal Ly/3 and 
N I A1134 contamination is marked in the FUSE spectrum. 



from the systemic velocity, but rarely redshifted and cer- 
tainly not anywhere near this much. Contrived models 
of He II damped Lya absorbers can be constructed at 
high redshift to account for an appar ent line-center shif t 
due to a strong red absorbing wing (|Zheng et al.|[200l . 
but no such model is possible with the low helium optical 
depth in the HS1700 sightline. It is even more puzzling 
that the FUSE spectrum, discussed below, shows no evi- 
dence for this excess, despite otherwise agreeing well with 
the COS spectrum. 

Therefore, we see no clear evidence for He II Lya emis- 
sion in HS1700. This leaves no detections of this line at 
good S /N, and only a few at low S/N and very low reso- 
lution (jSvphers et al.l [2009b). The lack of such emission 
despite predictions is puzzling, and may suggest that the 
models are overly simplistic. Because He II A304 scat- 
ters resonantly within the breadline cloud, some of the 
emission will be absorbed in the H I continuum, leav- 
ing only the surface layers to contribute to the observed 
line emission. Also, the Cloudy photoionization mod- 
eling code does not properly handle Bowen resonance- 
fiuorescence, where a He II Lya photon excites O III 
via a coincidence of transitions a round 303.8 A (jBowenl 
[19351: lOsterbrock fc Ferlalldl '2006. pp. 99-lOLl. How- 
ever, Cloudy does destroy these photons based on the 
presence of O III, which should suffice for our predic- 
tions. It is possible that the emission seen in HS1700 
could come from metals. Many models predict the blend 
of O III lines around A306 A to be the strongest metal 
emission anywhere near the break, although it is always 
predicted to be much weaker than He II Lya. The num- 
ber of quasars with good-quality spectra of the region 
near quasar-frame He II Lya continues to increase, and 
will hopefully offer more clues about both helium and 
metal lines in the EUV. 
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4.2. Comparison to FUSE 

HS1700 was obse rved by FUSE severa l t imes, most ex- 
tensive ly in 2003 (iReimers et al.ll2006( ). iFechner et all 
(|2006bf) performed the most detailed study of the FUSE 
data relating to the He II Lya forest. We present in 
Figure |3] an improved reduction of the FUSE spectrum 
(not previously published). The data were reduced with 
CALFUSE 3.0.7 using night-only data (288 ks) and PHA 
channels 4-16, with customized background subtraction. 
Individual exposures and segments were separately zero- 
corrected before coaddition. The region 1081-1087 A IS 
known to be unusually noisy, bu t we include i t in our 
figure for comparison with COS. IFechner et all (|2006bl ) 
excised this region from their analysis, but they used an 
earlier FUSE reduction with worse background subtrac- 
tion. 

Despite using night-only data, we still find noticeable 
geocoronal contamination in the Ff/6'i? spectrum, in part 
because of the large 30" x 30" (LWRS) aperture used. 
This is nearly 200 times the area (and hence sky back- 
ground) of the COS aperture, although both admit about 
the same amount of source flux. 

FUSE had R -20,000, much higher than COS G140L 
(i? ^ 2500), although it is possible that the resolution 
was degraded by pointing issues since the large aper- 
ture was used. As this is ideally high enough to re- 
solve the Lya forest, efforts were made to do so with the 
only other He II quasar observed w i th FUSE, HE234 7 
(|Kriss et al.l[200lt lzheng"erall l200l IShull et al.ll200l . 
Although the resolution may be good enough to do this 
rehably, the S/N likely is n ot, for either HE2347 or 
HS1700 f Fechner et aIll2006bD . Unfortunately our COS 
data cannot resolve the forest, although a newly avail- 
able bluer mode of the COS medium-resolution grating 
will allow such observations of some He II quasars down 
to z ~ 2.5. 

4.3. Comparison to STIS and the Intrinsic Spectral 

Index 

HS1700 was observed with STIS in 2003, for 2245 s 
using G140L {R - 1000) and 2166 s using G230L (i? - 
500). We plot this STIS data with our COS data super- 
imposed in Figure m HS1700 is know n to be high l y var i- 
able in the observed-frame FUV (jReimers et al l 120051 ). 
but we see only small differences in flirx and spectral in- 
dex between the two observations. Fitting a single power 
law below any strong pLLS Lyman breaks (A < 1900 A), 
we find that the STIS data give = 0.086+|^:^^ (for 

oc i/"-), while the COS data give = -O.lSt^Q H. 
However, because HS1700 is one of the few quasars 
with continuous spectral observations from the He II 
Lya break to the H I Lyman limit, it is possible to esti- 
mate all hydrogen absorption systems along the sightline. 
We performed a spectrum fit to the STIS data includ- 
ing pLLS, using the redshifts and column densities of 
IFechner et al.l ()2006af ) as initial param eters (many of the 
column densities are ultimately from IVogel fc ReimersI 
[19951 ). While the redshifts are derived from hydrogen 
and metal lines, and are therefore secure, the hydrogen 
column densities iVni are less certain. The low S/N and 
very low resolution (i? ~ 500) of the STIS NUV data 
makes fits problematic. Resolution is crucial even for 
measurements of the Lyman break, because the continua 



need to be well defined in a region with large amounts of 
contaminating absorption. 

The reasonably high S/N in the data allows us to 
detect pLLS absorption edges down to column densi- 
ties logA'ni > 15.2, corresponding to optical depths 
rm = (0.01)(7Vhi/10^^-2 cni-^). Our fit to the STIS 
data uses a single underlying power-law with 17 pLLS 
with logiVni > 15.2, nine of which have logiVni > 16.2 
(see Figure [SJ . For the mo st part we find colum n den- 
sities fairly close to those of IFechner et al.l (|200 6af). with 
the following exceptions. F or the z = 1.8450 syst em, we 
use logiVm = 16.80. The IFechner et al.l (|2006aD value 
of 16.21 was derived from fitting Lyman series lines, but 
is inconsistent with the strength of the break. For the 
z = 2.1989 system, we adopt l og TVht = 16.0 - while re- 
ported at logiVni = 15.44 in IFechner et al.l (|2006al ). a 
nontrivial break is helpful in the fit, although its column 
density is not well constrained. At z = 2.2895, we see 
no Lyma n break seen a ssociat ed with this system, de- 
spite the IFechner et al.l (|2006aD value of log A^m = 16.00 
(rm = 0.063). 

With the redshifts and column densities of the pLLS 
fixed, we then fit the COS data from the He II break to 
the G140L instrumental cutoff at 2000 A. We find the un- 
derlying EUV spectral indices to be a^, = —1.70 (COS) 
and ai, = —1.62 (STIS). These values agree with each 
other, considering systematic uncertainties. Quasars 
have widely varying spectral indices, but we note that 
this value is consistent with the average seen in radio- 
quiet HST EVY quasar spectra (a ^ = -1.5 7 ± 0.17, 
iTelfer et al.ll200l a^ = -1.41 ±0.21,lShull_eL.al. 2012b), 
though not with the FUSE average {a^ = -0.56t^J^, 
iScott et al.ll2004l ). The dramatic difference between the 
underlying spectral index and the apparent a,^ spec- 
tral index from a power-law fit to the COS data indicates 
the difficulty of discerning the true value without contin- 
uous spectral coverage from the Lyman limit. HS1700 
is, however, an extreme case with the number of pLLS it 
s hows. 

iGiroux et al.l ()1994[ ) analyzed the metal-ion absorbers 
in these pLLS, showing that their abundances were in- 
consistent with photoionization equilibrium. Instead, 
they proposed a multiphase medium with hot, collision- 
ally ionized gas producing much of the high ionization 
states of the C, N, and O ions. Howev er, the H I/He I 
ratios crucial to the IGiroux et al.l (|1994f ) analysis depend 
on FOS observations of He I that are not enti rely re- 
liable, as we note above. IFechner et al.l ()2006aD . using 
better data, fit the systems as photoionized only, but do 
not discuss the possibility of collisional ionization. 

We do not necessarily regard the fit of HS1700 with 
pLLS as superior to the simple power-law fit for calculat- 
ing Gunn-Peterson optical depths. The two continuum 
extrapolations are compared in Figure [51 The pLLS fit 
does account for the non-power-law nature of the con- 
tinuum, and should therefore indeed be superior, but 
there is nontrivial uncertainty in the column densities 
of the pLLS absorbers, and difficulty discerning the true 
continuum in data with low resolution and modest S/N 
(where th e true continuum i s easily confused with emis- 
sion lines: IShuU et al.ll2012bD . The high-transmission re- 
gion near 1040 A matches the power-law continuum quite 
well (Fig. [3]), and the pLLS fit is slightly unphysical there. 
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Figure 4. The 2011 COS (green) and 2003 STIS (black) spec- 
tra of HS1700, overplotted, v^ith error vectors plotted as dotted 
lines. The COS/G104L data are binned to one resolution element 
(0.56 A), while the STIS data are plotted at half a resolution ele- 
ment (0.58 A G140L, 1.54 A G230L). Both are corrected for Galac- 
tic reddening of E{B — V) = 0.026 using the extinction curve of 
ll'itzpatrick (1999 ). The close agreement in flux between the two 
observations is somewhat surprising, given the historically highly 
variable nature of this source. 




1500 2000 ^ 2500 

wavelength [A] 



3000 



Figure 5. The full STIS spectrum of HS1700 with an overplotted 
fit, including H I pLLS. For clarity, only the Lyman breaks are 
shown (tick marks for those with logA^jji > 15.4). Higher Lyman- 
series lines were included for the fit, improving the alignment near 
Lyman breaks, where the apparent break is slightly redward of 
the Lyman limit because of line overlap. Data from G140L (A < 
1715 A) are binned to one resolution element (~1.2 A), and data 
from G230L are binned to half a resolution element (~1.6 A). 



lying below the observed flux. The noise is fairly high in 
this region, so this is not a definitive test, but it does 
show that the pLLS fit is not necessarily better. The 
moderate-strength pLLS at z < 1 can be included or 
not — the COS data do not strongly show them or rule 
them out. Wiggles in the apparent continuum can be due 
to emission lines rather than rises and falls associated 
with pLLS. In this case, because we have metal absorp- 
tion data for these complexes (^Fechncr ct al. 2006a), it 
makes it more likely that they are truly pLLS. We do not 
see strong evidence for a s pectral index break at 2000 A 
(observed) , as suggested in iFechner et alj (|2006bf ) . 



Figure 6. A comparison of the pure power-law continuum fits to 
the best fit when including all the pLLS (though all but one of 
the Lyman breaks occurs at longer wavelengths than those shown 
here). The best-fit power-law is the solid black curve. The highest 
and lowest power-law fits for the 1000-1140 A region are dashed 
black lines (these are not the extremal fits at longer wavelengths). 
The best fit including all 17 of the logA'jji > 15.2 absorbers is the 
dot-dash blue curve. As one might expect, the fits are similar in 
the region where we have data, and differ noticeably only below the 
He II Lya break (1140 A). However, even at 1000 A, the difference 
is not significant for calculating Gunn-Peterson optical depths; see 
text for details. 



5. THE IGM IONIZING BACKGROUND 

One of the more interesting pieces of information in a 
spectrum of He II Lya, when used in conjunction with 
an H I Lya spectrum covering the same redshift, is the 
hardness of the EUV background. The ratio is charac- 
terized by 77 = A^Hcii/-^Hi — 4 X THeii/THi- While metal- 
line ratios from different ions can yield more detailed 
information, using helium gives the great advantage of 
almost complete coverage at every observable redshift. 
Very high values of rj (>500) are a sign of a soft ionizing 
field, dominated by stars, while low values suggest that 
the harder spectra of quasars dominates. 

Our plot of ri(z) is shown in Fig. [71 but unfortunately 
low resolution and the wavelength calibration uncer- 
tainty do not allow us to draw strong conclusions about 
specific features. The H I data are convolved with the 
COS LSF prior to calculating TcS, although this has lit- 
tle impact because each redshift bin contains at least 2.7 
COS resolution elements. 

We find a median (?7)incd — 31.7tg;4 (68% confidence 
interval from bootstrapping), with the 25*''-75*'' per- 
centile range being rj = 12-128. For the mean we find 
mean = 73t^°. IFechner et al.l ()2006b[ ) found {r])mcd = 
51 in the HS1700 FUSE data, and (77)niod = 79 when they 
restrict to a subset of the data where they argue thi is 
more robustly measured. Along the HE2347 sightline, rj 
at z = 2.4-2.73 was found to be (77)mod = 33 ( Shull et al 
I2OIOD . The mean was (77) = 78 ± 7 in iKriss et al. 



('20011 ) and simil ar in other works (jZheng et al.l I20(M 
Shull et al.ll201Cl although the latter quotes extremely 
large uncertainties). 

Most of the z = 2.2-2.7 IGM shows rj values consis- 
tent with quasars being the dominant source of ionizing 
photons, although some regions show softer backgrounds 
(77 > 500) that require strong filtering or local dominance 
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by softer sources such as star-forming galaxies. Such a 
hard background on average is expected in the He II post- 
reionization epoch, where there is no longer strong ab- 
sorption at the 4 ryd heli u m edg e, as there was at higher 
redshift. IMuzahid et al.l (|2011D claim that the average 
IGM region has a fairly high r] based on the HE2347 
sightline, and contend that low values {r] < 40) may 
be signs of gas that is not photoionized. The contribu- 
tion of coUisfonall y ionized He II was seen in simulations 
(jShuh et alJlMol ). Nonetheless, quasars are fully capa- 
ble of producing lower rj values, and indeed one might 
expect 77 30 for the average quas ar spectrum filtere d 
through an optically thin medium (jFardal et al.l 1 19981 ). 
aside from the many quasars that have spectra harder 
than average. In addition, median values of rj found 
here and some other works (e.g., ShuU et al. 2010) do not 
support the claim of high 77 being common after helium 
reionization, although we caution that we are using some - 
what lower resolution data than IMuzahid et al.l ( 2011t ). 
With Tqs this makes a difference, as higher-transmission 
regions dominate, but as this affects both hydrogen and 
helium spectra, it is unclear which way this would skew rj. 
Data in the He II Lya forest at good resolution and good 
S/N is the only way to conclusively examine this, but as 
yet no such data exist. By considering the "sawtooth- 
modulate d spectrum" prod u ced b y higher Lyman lines 
of He II, iMadau fc Haardtl (|2009 h predict a He II/H I 
ratio of 77 « 35 for optically thin filtering. 

Because 77 fluct uates strongly on fairly short scales 
(jShull et al.ll2010D . its measurement is sensitive to sys- 
tematic velocity offsets between the UV and the optical 
data. While we can confirm these averages for the post- 
He II reionization epoch, unfortunately, due to a larger- 
than-expected uncertainty in the wavelength calibration 
of the COS data, we are unable to speak about details 
in the 7/ fluctuations along the HS1700 sightline. 

Last we note a small but nonzero chance that there 
is a systematically incorrect estimation of our optical 
depths due to the possibility of wavelength calibration 
uncertainty in the sensitivity curve. In our determina- 
tion of the optical depths, we use the sensitivity curve 
of segment B data to predict the number of counts ex- 
pected for an unabsorbed spectrum. For A > 1080 A, 
the effective area of the COS detector changes rapidly 
(jMcCandliss et aLll2010D , and a small shift in wavelength 
could have an impact. For example, a wavelength shift 
of 2 A leads to an optical depth shift of At ~ 0.1. This 
shift, while not entirely trivial, is small compared to the 
values and other sources of uncertainty, so even in the un- 
likely case that such a large wavelength shift were present 
in the calibration data, its impact on our science would 
be minimal. 

6. PROXIMITY EFFECTS 

6.1. Line-of-Sight Proximity Effect 

Uncertainty in the systemic redshift of HS1700 has 
been a major source of uncertainty in interpreting the 
line-of-sight proximity effect. While nontrivial uncer- 
tainty remains, the redshift is now determined well 
enough to address this issue. Unfortunately, the COS 
data are of low resolution, and the break appears to be 
abrupt. We therefore consider the COS data in con- 
junction with the FUSE and Keck H I data, in Fig- 



ure [SI Wavelength calibration uncertainty (section 12. ip 
is larger than systemic redshift uncertainty for evaluat- 
ing the proximity effect, but the suggested zone is large 
enough — and well enough aligned with FUSE and Keck 
features — that this does not affect our results. 

At flrst glance, the FUSE data suggest there is essen- 
tially no proximity effect at all. This would be unex- 
pected for such a luminous quasar, although not unprece- 
dented as no proximity zone is seen in HE2347 either 
(but this is complicated by strong intrinsic absorption; 
iFechner et al. 2004; Shull et al. 2010). However, the ap- 
parent break in the FUSE spectrum is associated with 
a fairly strong absorber seen in H I as well. It is quite 
plausible that there is a He III zone surrounding this 
quasar out as far as the very high transmission in hydro- 
gen is seen (~14-24 comoving Mpc), interrupted by sur- 
viving dense absorbers. The COS spectrum, which does 
not suffer from geocoronal N I A1134 contamination like 
the FUSE spectrum does, supports this possibility. If 
the transmission peak near z ~ 2.73 is indeed associated 
with the proximity zone, then the highest-redshift point 
in Fig. [5] should be shifted up slightly, because some of 
this transmission lies in that bin. 

6.2. Transverse Proximity Quasars 

The transmission observed in the He II Lya forest is in 
part a reflection of the overall state of helium in the IGM 
at a given redshift, a measurement of the global ioniz- 
ing background that evolves during helium reionization. 
However, it also can depend on local ionizing sources. In 
particular, a quasar that is near the sightline can affect 
the ionization state of the gas along the sightline, lead- 
ing to a transverse proximity effect. Measurements of 
this effect can give information about quasars, notably 
lifetimes and b eaming angles (Furlanctto & Lidz 201l|; 
ILu fc Yull20ll . 

Measurement of increased transmission due to the 
transverse proximity ef fect on the H I Lya for est has 
been unsuccessful (e.g.. iKirkman fc Tvtleii |2"008[ ) . likely 
due to factors such as quasars existing in regions that 
are denser than average, counteracting their increased 
ionizing flux. It has proved much more fruitful to look 
at ion ratios, which give information on spectral hard- 
ness and break the degeneracy between ionizing flux and 
IGM density. Metal syste ms can be used for such ra- 
tios ('Gon calves et al.ll2008f) , but the vast majority of the 
path length through the Lya forest has no observable 
metal systems. It has therefore been rewarding to com- 
pare the hydrogen and helium Lya forests, which gives 
a measurement of spectral hardness that can be used 
anywhere. Several detections of transverse proximity 
quasars have been made thi s way, for the sightl ines to 
00302-003 (Jakobsen et al. 2003; Worseck fc Wisotzkil 
12006.') and HE2347 (Worseck ct a l. 200711 " 

We present the first survey for high-redshift quasars 
near the HS1700 sightline (Figure O Table [I]). The field 
near HS1700 is covered by SDSS imaging, with good- 
quality ugriz data. Unfortunately, the region nea r z ^ 
2.7 is a known redshift desert for SDSS (Richar ds et al.l 
I2009f ). due to quasars overlapping the stellar locus in 
color-color space. Efforts have been made to construct 
photometric quasar catalogs none theless, and in our se - 
lecti on we draw on tw o of these, iRichards et al.l (|2009[ ) 
and iBovv et all (|2011[ ). Because the targets observed 
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Figure 7. Optical deptlis of He II and H I Lya as a function of redshift, with 4 X 'rHcii/''"Hl — V plotted in the lower three planels. The first 
Tj panel uses the wavelengths as reported by CALCOS, while the second shifts by A = 2 A and the third by A = 3.4 A, reflecting uncertainty 
in the wavelength calibration. Our best estimate of the true shift is 2.0 A; see text for details. The data are binned to Az = 0.005 (about 
2.7 COS resolution elements) for z > 2.4, and Az = 0.01 for z < 2.4 due to lower S/N in this region. Errors shown are 68% including 
Poisson uncertainty in FUV counts and continuum uncertainty. Uncertainty in thi eff is smaller than the plotting symbols. The optical 
depth adjustments from accounting for contaminating (non-helium) absorption in the FUV spectrum are included, but the original values 
all lie within the plotted confidence intervals. 
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Figure 8. Normalized flux for the COS, FUSE, and Keck spectra, shifted to redshift space for He II and H I Lyo. The redshift of HS1700 
is marked with the dashed vertical line, with approximate redshift uncertainty shown in the COS panel. Error vectors are overplotted 
as dashed lines for the UV data. COS data are binned to ~85 km s~^ (just over half a resolution element), FUSE data are binned 
to ~27 km s~^ (1^2 resolution elements), and Keck data are binned to ~7 km (one resolution element). The FUSE spectrum is 
contaminated by geocoronal N I A1134, which is marked. 
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spectroscopically were chosen for a variety of factors, we 
emphasize that this survey has considered most "good" 
candidates within ~20' (~10 proper Mpc) of the sighthne 
with r < 22, but it is not complete. 

We note that more efficient differentiation of quasars 
and stars can be achieved with m any epochs of imagin g 
allowing variability selection fe.g. JMacLeod et al.ll201l"l ). 
or with IR bands for those quasar s bright enough to be 
detected in WISE ^Wu et al.ll2012D . Slitless spectroscopy 
is difficult for large fields in the optical, due to crowding, 
but it has been successfully used to iden tify transverse 
proxi mity quasars in some southern fields (jWorseck et alJ 
[2001 . 

Our spectra of candidate transverse proximity quasars 
are shown in Figure |9l To verify candidate quasars and 
obtain redshifts, we used the Dual Imaging Spectrograph 
(DIS) on the 3.5m telescope at Apache Point Observa- 
tory (APO). DIS has a blue side and a red side, with 
a constant dispersion of 1.8 A pixel^^ in the blue and 
2.3 A pixel" ^ in the red, with two pixels per resolution 
element. We use the blue side for A < 5425 A, and 
the red side longward of this. Independently calibrated, 
the two sides normally match fluxes seamlessly, but fac- 
tors such as a poor trace of a faint spectrum can cause 
matching problems such as the one seen in Figure [HI for 
SDSSJ170219+640454. 

The identification of SDSSJ170032-t-640524 is tenta- 
tive. There is only one very clear emission line in this 
spectrum, and the redshift depends on identifying this as 
Lya. C IV 1549 does appear to be present at a consis- 
tent redshift, but no other lines are seen. If the strong 
emission line is instead Mg II, we should see II/3 in the 
spectrum, but we do not. Spikes near 5600 A and 6300 A 
are narrow — consistent with purely instrumental broad- 
ening (~250 km s~^), and thus likely noise — while the 
identified emission line has a FWHM of '^1300 km s~^. 
This would make it narrow for a quasar, but much too 
broad to be from a star-forming galaxy. 

SDSSJ170335+642603, at z = 2.381, is unfortunately 
of minimal observational interest for transverse proxim- 
ity effects. At that redshift. He II Lya is observed at 
1027 A — right on top of geocoronal H I Ly/3, a strong 
contaminating feature in our COS spectrum. 

There is a significant protocluster of galaxies at z ~ 
2.30 within a few arcminutes of the sightline; this 
protocluster should generate a soft UV background, 
with many star-forming galaxies but no reported AGN 
([Stcidcl et al. 2005). This could therefore potentially 
affect H I Lya transmission in the sightline, but not 
He II Lya transmission. However, we see no evidence 
for a soft photoionization effect (Fig. [7]). 

7. CONCLUSIONS 

We have presented a new COS spectrum of the bright 
He II quasar HS1700. This spectrum has the highest S/N 
and resolution of any UV spectrum of this quasar in the 
region 1170 < A < 2000 A, and higher S/N and much 
better background subtraction than the FUSE spectrum 
at shorter wavelengths. This allows us to accurately 
calculate Tcff(z) for this sightline, finding values signif- 
icantly higher than previously reported, at in agreement 
with values found for the HE2347 sightline at the same 
redshifts. The He II post-reionization epoch is charac- 



terized by steadily falling effective optical depths, and a 
UV background that is varying but consistent with being 
dominated by the hard spectrum of quasars. 

HS1700 was one of the few He II quasars with a claimed 
detection of He II Lya emission, and by far the bright- 
est. Our COS spectrum does indeed show an excess of 
flux near the He II Lya break, significantly above any 
reasonable continuum fit. However, the apparent peak 
of this flux is about 3000 km s~^ away from He II Lya at 
the quasar rest-frame, even when using the much higher 
and more accurate systemic redshift of zhsitoo — 2.75. 

Transverse proximity quasars are difficult to detect, 
but one of the most promising methods has been to look 
at variations of the spectral hardness using rj. We present 
four candidate proximity quasars along the HS1700 sight- 
line. Although none are clearly detected in the present 
data by r/ variations, two are at high enough redshift 
that they could be observed at i? > 10,000 with a newly 
available COS mode. This resolution, combined with 
the much greater wavelength calibration certainty com- 
pared to the present G140L data, could definitively an- 
swer whether or not they are observed in the transverse 
proximity effect. Either answer has interesting implica- 
tions for quasar physics. 

Background determination is one of the largest sources 
of uncertainty for studies of effective optical depth. In 
this paper we present and use for the first time a new 
method of COS background subtraction that substan- 
tially improves on the pipeline strategy. This method 
will be even more important when looking both at fainter 
targets, and crucially, prior to the end of helium reion- 
ization, when the optical depths are large. 

We thank Rob Simcoe for providing the Keck spec- 
trum of HS1700, Jerry Kriss and Jennifer Scott for pro- 
viding their final reduction of the FUSE spectrum, and 
Justin Ely for useful conversations regarding the COS 
background. 
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NNX08AC146 and NAS5-98043 to the University 
of Colorado at Boulder. 
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Table 1 

APO/DIS Observations of Potential Transverse Proximity Quasars 



Name 




R.A. 


Decl. 


Redshift 


Impact Parameter 




Exp. Time 


Obs. Date 






(J2000) 


(J2000) 




[proper Mpc] 


[mag] 


[min] 


(UT) 


SDSSJ170032.86- 


h640524.8 


255.13694 


64.09025 


2.277 ± 0.012*' 


3.72 


21.7 


60 


2011.08.27 


SDSSJ170335.96- 


h642603.8 


255.89987 


64.43439 


2.381 ±0.003 


10.9 


21.7 


90 


2011.05.31 


SDSSJ170029.48- 


[-635750.9 


255.12286 


63.96415 


2.588 ±0.014 


7.40 


21.2 


60 


2012.05.29 


SDSSJ170219.00- 


h640454.5 


255.57921 


64.08182 


2.625 ±0.005 


5.64 


20.9 


100 


2012.05.16 



From SDSS DR8 IIEisenstein et al.|[20Tl)). These are not n ecessarily the magnitudes used in constructing the photometric 
catalogs — for example, in DR7 (used by [Richards et al."2009) the z = 2.381 quasar was 0.8 magnitudes brighter in r, due to 
some combination of quasar variability and altered data processing. 
^ Line identifications are tentative; see text for details. 
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APPENDIX 

DETERMINING THE COS FUV BACKGROUND 

The background of COS is dominated by dark current from the detecto r, except at specific wavelengths where geo- 
coronal hne emission dominates. Scattered and Zodiacal light is neglig ible (|Svphers et al.llMl . The current CALCOS 
(2.18.5) method for background determination uses two windows offset from the primary science aperture (PSA) in the 
cross-dispersion {y) direction, with a runni ng average ove r 100 p ixel£| in the dispersion (x) direction. The drawbacks 
of this method are discussed extensively in lSvphers et al.l (|2012( ). and here we present the first implementation of the 
improved method suggested in that work. 

The primary problem is so-called y-dip, where counts in the PSA have lower pulse-height amplitudes (PHAs) than 
those in background regions (see Fig. 6 in Syphers ct al. 2012). Because of detector thresholds and PHA cuts, this 
means the background in the PSA is actually noticeably lower than the background inferred from regions offset in y. 
(We note that this problem is lessened but by no means removed if no PHA filtering is imposed, and such inclusive 
analysis adds a large amount of noise.) On July 23, 2012, the PSA was shifted to a new lifetime position on the detector, 
but y-dip will again become an issue as the new position is used more. Also, an underestimate of the background could 
o ccur if backgrq t ind w indows overlap regions previously affected by y-dip. 

iSvphers et all (|2012[ ) suggest using data from the COS FUV dark monitoring programs (GO 11895, 12423, and 12716) 
to precisely characterize the dark current in the PSA, and use this for background subtraction. Although ideally the 
dark would be recorded without binning, the very low dark current precludes this. The rate at 2 < PHA < 30 is 
~ 2 X 10~^ counts pixel"^, and the dark monitoring programs collect data fairly infrequently — once a week for 
about 27 ks per month for 11895 (September 2009 to October 2010) and once every two weeks for 13 ks per month for 
the more recent programs. At this pace it would require years of collection to obtain on average a single background 
count per pixel. The y-dip evolves noticeably with time, because it is caused by burn-in of the detector, so averaging 
data over too-long intervals is problematic. Absent more frequent dark monitoring, we must both average moderately 
long time intervals and smooth the background. 

We exclude dark exposures with unusually high overall count rates (we used the limit 6 x 10~^ counts s""'^ pixel"^, 
but caution that the average background does vary noticeably with the solar cycle). There are two categories of 
individual pixels with excessive counts: transient hot spots and consistent hot spots. Individual exposures are checked 
for transient hot pixels, and we exclude any pixel with an excursion that would be expected < 1% of the time on 
the detector, assuming Poisson noise. Consistent hot pixels also exist and are excluded. These show up as very 
improbable excursions on the coadded master dark, but on any individual exposure they are often probable enough 
to not be excluded. A true pixel-by-pixel master dark would include such elevated rates, but our smoothed version 
does not, as that would contaminate surrounding pixels. In addition, when working in regions that are very sensitive 
to background rates, we recommend excluding regions known to have elevated backgrounds (data quality flag 32). 

We also impose a PHA cut on our data that is more stringent than the pipeline standard. The CALCOS pipeline 
by default includes 2 < PHA < 30 (recorded PHA ranges from to 31). However, the PHA distribution for sources 
is sharply peaked at lower PHA, while the PHA distribution for dark counts is much flatter. For the segment A, tests 
with HS1700 showed that source counts contributed signiflcantly only to the bins 1 < PHA < 14, and there were 
no detectable source counts for PHA > 16 (whereas about 30% of dark counts have PHA > 16). The detection of 
source counts in the PHA = 15 bin is marginal, and in any case < 0.07% of source counts have this PHA. Counts with 
PHA = 1 are a more subtle case, as ~70% of them are source counts, but they are not included in the flux calibration. 
We choose to neglect these counts — only ~0.5% of source counts lie in this bin, and the important thing for our optical 
depth analysis is to predict source counts for an unabsorbed spectrum (using the flux calibration to convert from 



^ We refer throughout this discussion to pixels, which unless 
otherwise specified means the two-dimensional pixels created by 



digitizing the output of the cross delay line detector. While the 
COS detector itself does not have physical pixels, they are funda- 
mental in the data the observer receives. 
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extrapolated flux to counts) and measure them in exactly the same way. Segment B is similar, although at slightly 
higher PHA due to the higher voltage of that segment. There are no detectable source counts with PHA = 2 3, 
but nontrivial numbers have PHA = 15-16. (We note that there are substantial source counts with lower PHA in 
those portions of the detector strongly aflFected by gain sag due to geocoronal Lya, but all such regions lie below the 
Galactic cutoff of 912 A for G140L data.) Our analysis of the HS1700 data is therefore restricted to 2 < PHA < 14 on 
segment A, and 4 < PHA < 16 on segment B. We do recommend analyzing a specific data set prior to making PHA 
cuts, because the PHA distribution has several time dependencies, notably including loss of sensitivity in the exposed 
portion of the detector, shifts to a new lifetime position (as was done in July 2012), or changes in the high-voltage 
settings. 

The spectrum extraction width (in the cross-dispersion direction) also contributes substantially to how much back- 
ground is included. The default extraction width for COS G140L data is 57 pixels in y, but the spectrum width is 
much smaller than this, and therefore many unnecessary dark counts are included. For the HS1700 data analyzed here, 
we use an extraction width of 30 pixels for segment A, and 35 pixels for segment B. This is fairly generous, as >95% 
of all source counts are contained within 20 pixels over most of the detector, but we err on the side of being inclusive 
because the flux calibration will be off if too many source counts are cut out. Also, the spectrum widens slightly at 
detector edges, including high wavelengths on the segment A (A > 1900 A) and the usable portion of the segment B 
(1000 < A < 1175 A). With our chosen widths, we preserve essentially all source counts for most of the spectrum, and 
>98% at all wavelengths. 

We directly verify fiux accuracy by comparing standard extracted spectra to our extracted versions, and finding 
essentially identical flux levels in regions where the background is negligible. With both PHA cuts and a narrower 
extraction width, we reduce our background to about 40% of what it would be in a default extraction, without 
appreciably losing source counts. 



